Search for: All records

One of the most fundamental problems that has no efficient solutions on classical computers is simulation of quantum systems. It has been long hypothesized that quantum computing devices are naturally more suitable for this task, but many aspects of practical implementations of such simulations remain unknown. One particularly important kind of these simulations is the simulation of molecular dynamics, i.e. prediction of time evolution for a system of interacting particles. In this work we show how a quantum annealer can be used to carry out such simulations by solving differential equations of motion, on the example of the hydrogen molecule. Although the considered system is simple, our method is well scalable and can be readily applied to more complicated systems as annealers with larger number of qubits become available. Importantly, the method is general and can be used to solve arbitrary systems of ordinary non-linear differential equations, which can be helpful not only in the field of computational chemistry, but in many other fields as well. 

A theoretical approach is developed for the description of all possible recombination pathways in the ozone forming reaction, without neglecting any process a priori, and without decoupling the individual pathways one from another. These pathways become physically distinct when a rare isotope of oxygen is introduced, such as 18O, which represents a sensitive probe of the ozone forming reaction. Each isotopologue of O3 contains two types of physically distinct entrance channels and two types of physically distinct product wells, creating four recombination pathways. Calculations are done for singly and doubly substituted isotopologues of ozone, eight rate coefficients total. Two pathways for the formation of asymmetric ozone isotopomer exhibit rather different rate coefficients, indicating large isotope effect driven by ΔZPE-difference. Rate coefficient for the formation of symmetric isotopomer of ozone (third pathway) is found to be in between of those two, while the rate of insertion pathway is smaller by two orders of magnitude. These trends are in good agreement with experiments, for both singly and doubly substituted ozone. The total formation rates for asymmetric isotopomers are found to be somewhat larger than those for symmetric isotopomers, but not as much as in the experiment. Overall, the distribution of lifetimes is found to be very similar for the metastable states in symmetric and asymmetric ozone isotopomers. 

Scattering resonances above dissociation threshold are computed for four isotopically substituted ozone species: 16 O 18 O 16 O, 16 O 16 O 18 O, 18 O 16 O 18 O and 16 O 18 O 18 O, using a variational method with accurate treatment of the rotation–vibration coupling terms (Coriolis effect) for all values of the total angular momentum J from 0 to 4. To make these calculations numerically affordable, a new approach was developed which employs one vibrational basis set optimized for a typical rotational excitation ( J , Λ ), to run coupled rotation–vibration calculations at several desired values of J . In order to quantify the effect of Coriolis coupling, new data are contrasted with those computed using the symmetric-top rotor approximation, where the rotation–vibration coupling terms are neglected. It is found that, overall, the major properties of scattering resonances (such as their lifetimes, the number of these states, and their cumulative partition function Q ) are all influenced by the Coriolis effect and this influence grows as the angular momentum J is raised. However, it is found that the four isotopically substituted ozone molecules are affected roughly equally by the Coriolis coupling. When the ratio η of partition functions for asymmetric over symmetric ozone molecules is computed, the Coriolis effect largely cancels, and this cancelation seems to occur for all values of J . Therefore, it does not seem grounded to attribute any appreciable mass-independent symmetry-driven isotopic fractionation to the Coriolis coupling effect.